Abstract 



The general hydro-thermodynamic system of equations in 2+1 dimen- 
sions with arbitrary equations of state (Taylor series approximation) is 
split to eigen modes: Tollmienn-Schlichting (TS) wave and two acoustic 
ones. A mode definition is realized via local relation equations, extracted 
from the linearization of the general system over boundary layer flow. 
Each such connection defines the subspace and the corresponding pro- 
jector. So the division is performed locally and could help in initial (or 
boundary) problems formulation. The general nonlinearity account de- 
termines the specific form of interaction between acoustic and vortical 
boundary layer perturbation fields. After the next projecting to a sub- 
space of Orr-Sommerfeld equation solution for the TS wave and the cor- 
responding procedure for acoustics, the equations go to one-dimensional 
system that describes evolution along the basic stream. A new mechanism 
of nonlinear resonance excitation of the TS wave by sound is proposed 
and modelled via four-wave interaction. Subjectclass: Primary 35Q30 ; 
Secondary 76F70 , 35Q72; Keywords: keywords fluid mechanics, bound- 
ary layer, projector to eigen modes, Tollmien-Schlichting waves, acoustic 
waves, nonlinear resonance, N-wave system. 
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1 Introduction 

It is well-established fact now that a free-stream and a surface disturbances 
strongly affect the processes in boundary layer which compose a transition to a 
turbulent state. This transition in turn determine such important parameters 
of a rigid body in fluid mechanics as skin friction and heat transfer. The clas- 
sification of the free stream was done first by Kovasznay et al |jj . It is shown, 
that within linear approximation a general small-amplitude unsteady distur- 
bance in the free stream can be decomposed into three independent different 
types: acoustic (A), vortical and entropy modes. Only the first one relates to 
pressure fluctuations propagating with the sound speed, the last two don't case 
any pressure perturbation. The common idea of many investigations is to pick 
out length and time scales of each of these disturbances that would make the 
possibility of the Tollmienn-Schlichting (T-S) waves nonlinear generation 
This way the mechanism of T-S wave generation by convecting gusts interact- 
ing with sound was justified and numerically investigated H . Appropriate scales 
of boundary roughness are proved to generate T-S waves both theoretically and 
experimentally Q. 

In spite of abundant efforts ( see the big introduction and the citations in 
[||), devoted to the problem of a search of an effective control mechanism that 
support a cumulative direct energy exchange between T-S an acoustic modes. 
Let us mention three important papers about the general, local and distributed 
acoustic receptivity ||, ||, @- The results however do not look complete: 
we revisit the problem in the all-perturbations approach 0, |) starting from 
boundary layer (BL) as a background. The perturbations are considered only 
over the stationary boundary layer, we do not account here (but plan to do 
it) the layer field as a dynamic variable. This scheme gives a possibility to 
study mutual interactions on the base of model integrable nonlinear evolution 
equation. In our approach the description is local by the construction and do 
not need averaging procedure pOj. It covers the known results and give new 



hopes for understanding related phenomena appeared in papers from 11 to 
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[O . We would note that the initial stage of the process of its structure changes 
manifests. This effect is similar to heating or streaming generation by acoustic 



waves 13 , it is a development of the stationary mode and corresponds to an 
initial stage of the BL reconstruction. 

In this paper we concentrate our efforts on the mathematical formalism: in- 
troducing the complete set of basic modes we transform the fundamental system 
of standard conservation laws of fluid mechanics to a set of equivalent equations. 
In linear approximation the specific choice of new dependent variables split the 
system to the set of independent equations for the given modes, the account of 
nonlincarity naturally introduce the interaction by projecting the fundamental 
equations set in a vector form. Going to the nonlinear description, we use iter- 
ations inside the operator by the small parameter related to amplitude (Mach 
number for acoustics) and viscosity (Reynolds number). We also analyze the 
possibilities of resonant interaction of quasi-plane waves on the level of so-called 
N-wave systems § . Being integrable such systems admit explicit solutions and 
plenty of conservation laws. Hence the detailed investigation of situation is 
possible in this approximation. 

The mentioned types of waves (T-S and A) are defined by eigenvectors of 
the linearized system of dynamic conservation equations for the free flow. Once 
defined, eigenvectors (or modes) are fixed and independent on time. The process 
by which the free stream disturbances are internalized to generate boundary- 
layer instability waves is referred to as receptivity. The basic idea of this paper 
is to define the T-S and acoustics modes as eigenvectors of the same system 
for a viscous flow over a rigid boundary. The eigenvectors of the viscous flow 
go to the known limit in the free stream over a boundary. Our idea is to fix 
relations between specific perturbations (velocity components, density and pres- 
sure) for every mode followed from the linear equations, to construct projectors 
and applied them for the nonlinear dynamics investigation. Thus, in the linear 
dynamics the overall field may be separated by projectors to the independent 
modes at any moment. Fixing the relations when going to the nonlinear flow, 
one goes to a system of coupled evolution equations for the modes. These ideas 
have applied to nonlinear dynamics of exponentially stratified gas and bubbly 
liquid |l and other problems ||. 

The principal difference between this consideration and previous ones is the 
necessity of variable coefficients incorporation that arise from a boundary layer 
structure. It means that such coordinate dependence impact the projectors 
structure: the projector operator should be constructed by nonabelian entries. 
By other words the matrix elements of the projector matrix will be operator- 
valued ones. In fact we revisit and develop the first attempt in this direction 
that had been made recently |l4| , jl5| . 

The T-S wave takes the place of a vortical mode in unbounded space, the 
difference is due to the linearization on the different background. When a sta- 
tionary boundary-layer flow like Blausius one appears, a linearization should be 
correctly proceeded with account of the boundary-layer flow as a background, 
that would lead obviously to the other features of the vortical mode then that 
in unbounded space. An important feature of T-S mode is non-zero disturbance 
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of pressure already in the linear theory. In the three-dimensional flow, there 
exist two T-S modes, two acoustic ones (corrected by background flow), and 
one entropy mode as well. 

The equations of interaction are derived in Sec. 3 by means of the division 
of the perturbation field on these subspaces and projecting the system of the 
basic equations on the same subspaces. This transformation that in fact is 
nothing but a change of variables allows to proceed in a choice of adequate 
approximation. Moreover we could analyze possibilities of a direct nonlinear 
resonance interaction account. The results are the following: two- wave and 
three-wave interaction does not contribute due to the structure of the interaction 
terms in its minimal possible order. Hence (in this order) only the four-wave 
interaction display a resonance structure. We derive the correspondent four- 
wave equations and analyze it in Sec. 4. 

2 Basic equations treating the equations of state 
in the general form. 

The mass, momentum and energy conservation equations read: 

%+V{pv) = i) (1) 



- + (vV)v 



-Vp + i]Av +(•?+!) V( W) 
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Here, p,p are density and pressure, e, T - internal energy per unit mass 
and temperature, r?, <r, are shear, bulk viscosities, and \ ~ thermal conductivity 
coefficient respectively (all supposed to be constants), v is a velocity vector, 
Xi - space coordinates. Except of the dynamical equations (1,2,3), the two 
thermodynamic relations are necessary: e(p, p),T(p, p) . To treat a wide variety 
of substances, let us use the general form of the caloric equation (energy) as 
expansion in the Taylor series: 
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The background values for unperturbed medium are marked by zero, perturba- 
tions of pressure and density are denoted by the same characters (no confusion is 
possible since only perturbations appear below ), C v means the specific heat per 
unit mass at constant volume, E±, ...Oi, ... are dimensionless coefficients. The 
two-dimensional flow in the coordinates x (streamwise distance from the plate 
(model) leading edge), z (wall- normal distance from a model surface) relates to 
the two-component velocity vector 

v = (u,w) + u , (6) 

where Hq = (Uq(z)^O) denotes the background streamwise velocity and (u,w) 
stands for velocity perturbations. The system (|l| - ||) with account of (|], ||, || ) 
yields in 

(< + §)(& + fe) = ~ + £§f 

^(fr + ^oiO + fi-ryA^ + f) (fe 

-p^-p^+i-M 

where the constants Z and <S* are defined by 

Z= (-1 + 2^E 3 + E^ /E u 
S = ( 1 + E * + 2i? 4 + ^ E 5 ) 

The constant c = \J M p~ E ^ 2) has the sense of linear sound velocity in the medium 
under consideration when Uq = 0. The right-hand side of equations involves the 
quadratic nonlinear and viscous terms as well as linear ones related to ther- 
mal conductivity, no cross viscous-nonlinear terms accounted. In fact, the third 
equation follows from the energy balance and continuity equation. No assump- 
tions on flow compressibility was not done yet. 



d 2 w 



(7) 



OZ ' 
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3 Modes in the linear approximation. 

The basic system (Q-|^) contains four dynamic equations and therefore, there 
are four independent modes of the linear flow: two acoustic ones, vorticity and 
heat modes. This is a classification by Kovasznay who defined the acoustic 
modes as isentropic and irrotational flow, and refers to the two last ones as to 
frozen motions, the vorticity one relating to the absence of pressure and density 
perturbations and the heat one relating to the very density perturbation. The 
system (|l}||) involves three equations indeed, therefore only the three modes 
may be extracted - two acoustic and vorticity ones. That is due to the structure 
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of the heat mode where the only density perturbation occurs. Strictly speaking, 
the term treating a thermal conductivity in the third of equations (^) includes 
the density perturbations and the linearized system (^) (that defines the modes 
as possible types of flow) is not closed. If there was no this term at all there is a 
simple linear relation between density and pressure perturbations for the both 
acoustic modes. The presence of thermal conductivity corrects this relation as 
it was shown in p8| . When the effects of thermal conductivity are small, the 
corresponding terms may be placed to the right-hand side of equation together 
with nonlinear ones and be accounted further. So the excluding of the dynamic 
equation for density serve just a simplification of a problem suitable in the view 
of its extraordinary complexity when the heat mode is out of the area of interest. 

The main idea is to define modes accordingly to the specific relations of the 
basic perturbation variables following from the linearized system of dynamic 
equations (|l|-||). In general, the procedure is algorithmic and may be expressed 
as consequent steps : to find the dispersion relation and its roots that determine 
all possible modes, and to define relations between specific variables for every 
mode. Later, projectors follow immediately from these relations , they separate 
every mode from the overall field of the linear flow exactly and serve as a tool 
for the nonlinear dynamics investigation. 

Using the left-hand part of the system (^) as a basis of modes definition 
and introducing the two non-dimensional functions Vq(z) — Uq(z)/U oq , 4>{z) = 
Vo z (z)Iq , we rewrite it in the non-dimensional variables 

= x/Iq,w* = w/Uoo,u* = u/£/oo,t* = tUoo/l ,p* =p/p a U^ . (9) 

The value C/qq marks velocity of a flow far from the boundary, and Iq - boundary 
layer width. In the new variables (asterisks will be later omitted) ([?]) with zero 
right side reads : 

dp dp _ 2 ( ® u dw\ 
dt dx \dx dz J 

du Tr du , dp _ i * n-i f d 2 u d 2 w \ „ 

1 - + V — + cbw + 7 f-Re 1 Au-R 1 l— + ——)=0 11 
at ox ox \ ox dxdz J 

& w _^_y d w _|_ jj, ( ° 2 u ^ d 2 w\ q 

dt dx dz \dxdz dz 2 J 

with parameters e = Uqo/c (the Mach number), Re — Uoclopo/rj is the Reynolds 
number base on the length scale, and R = C/oo/oPo/ + ?) • 

3.1 The Tollmienn-Schlichting mode. 

Formally, the limit of incompressible fluid ( e = ) corresponds to the vorticity 
mode. The first relation for velocity components is well-known [pi: 

du/dx + dw/dz = 0. (13) 
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From the equations (]Tl],[T^) an expression for the pressure perturbation fol- 
lows: 

2(j)dw/dx + Ap = (14) 

where A = d 2 jdz 2 — k 2 stands for the Laplacian so far the d/dx equivalent 
operator (multiplier) — ik is used. Both ( |l3| ) and (|l4| ) define the TS mode due 
to relations of the specific perturbations of pressure and velocity components. 
Since the geometry of the viscous flow over boundary supposes strong non- 
uniformity in the vertical direction, all disturbances may be thought in the 
basis of not plane waves but in the functions like ^(x, z) = ip(z) exp(iojt — ikx) 
. So, one has to leave vertical derivatives that usually are large in comparison 
with the horizontal ones. The result is rather obvious hence we do not introduce 
a special small parameter. The vector of the T-S mode may be chosen as: 



Pts (15) 



Note also, that the equations ( |10| - |12| ) yield the well-known equation for the 
TS mode JT^] , when rewritten for the new variable such the stream function 
( u = d^/dz , to = —d^/dx, with the obvious restriction to the solenoidal 
velocity field (|l|): 

Ad^/dt + VoAdV/dx - d^/dx ■ 84>/dz - .Re^A 2 * = (16) 

The well-known Orr-Sommerfeld (OS) equation follows from (^6|): 





( PTS \ 
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UTS 


-1 




V W TS J 






[Vq(z) - c] [d 2 ip/dz 2 - k 2 ip] - ^d<j)/dz = — - [d^j/dz^ - 2k 2 d 2 ^/dz 2 + k^] 

(17) 



That equation is an initial point of the laminar flow stability theory and for 
every pair (k,Re) determines an eigenfunction ip( z ) an( i complex phase velocity 
c = uj/k — c r + ici . The sign of Cj is namely a criterion of flow stability: a 
negative value corresponds to the growth of perturbation and therefore to the 
non-stability of the flow. 

3.2 Acoustic modes 

The potential flow imposed two acoustic modes with 

du/dz — dw/dx = 0. 

In the limit of Re^ 1 = , Br 1 =0,0 = 0, (|To|- 12) naturally goes to the acoustic 
modes dispersion relation that is directly connected with the wave operator. 
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We would not consider here the perturbation velocity field changes, forced by the 
ambient movement of the fluid. It could be account by the perturbation theory 
to be developed here. Then, two acoustic modes are defined with relations 
between specific perturbations: 
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(19) 

Here the square root of the operator A is defined as integral operator via 
Fourier transform. 



4 Projectors in the linear problem. 

The TS and two acoustic modes are determined by relations of specific perturba- 
tions or, in the vector form ([l5|), ([l9|). The superposition of such disturbances 
appear as all possible types of flow except of the heat mode. Every mode is 
completely defined by one of specific perturbations - pressure or velocity com- 
ponents since there are strict and local relations between them. Practically, in 
a linear flow, the overall perturbation may be de-coupled into modes by the 
corresponding orthogonal projectors. The arbitrary perturbation then is a sum 
of modes which, taking into account (R~5|) , (ft9|) looks: 
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H 







PAl + PA2 + PTS \ I PA1+ PA2 + PTS 

UA1 + U A 2 + u T s = HpAl - Hp A 2 + Kp TS j (20) 
WA1 + WA1 + WTS ) \ MpAl - Mp A 2 + QPTS 



with operators 



M = eA-V 2 d/dz, 

K = -^d/dz^A, W 

Q = ih A 

The link ( po|) may be considered as a one-to-one map of dynamical variables 
that immediately yields in the case of TS wave the projector : 

/ -2fc 2 A- 1 0'9/9zA- 1 -2ifc 3 A- 1 0A" 1 

P TS = d 2 /dz 2 A~ 1 ikd/dzA- 1 j (22) 

\ ikd/dzA- 1 -k 2 A~ l . 

For a right and left acoustic waves one has 

Pax = 
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(24) 



The projectors possess all properties of orthogonal projectors and their sum 
is unit matrix since all eigenvectors of linear system are accounted. 



p2 



P TS , Pts + Pai + Pa2 = I, Pts ■ Pai = Pts ■ Pa2 



0, 



where / and are unit and zero matrices correspondingly. We would stress 
that the operators defined by the ( ]22] , jgjj , 0) contain operator-valued matrix 
elements. Hence the derivation of their explicit form as well as the properties 
mentioned above could be checked taking the nonabelian nature of the operators 
into account. In the linear flow, the projectors separate every mode from the 
overall perturbation, for example: 
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(25) 



and so on. Moreover, acting by a projector on the basic system of dynamic 
equations (|l0|-p"2|) , yields a linear evolution equation for the mode assigned to this 
projector. One produces three equations indeed for every specific perturbations 
which are essentially the same with account of relations (|l4|) , ( |l9| ) . For the first 
rightwards progressive acoustic mode the evolution equation reads: 



dpAi/dt + V a d P Ai/dx + e-^V = 



(26) 



The equation for the second (opposite directed) acoustic mode is produced by 
acting the projector pA2 on the basic system and differ from ^6] only by the sign 
before the last term. Combining equations for these directed acoustic modes 
one arrives at the wave equation of second order 

(d/dt + Vod/dxfp - e~ 2 A P = 0. 

This equation appears as a limit of more general one: that relates to both 
acoustic modes and can be found in (lGl. 
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5 Nonlinear flow: coupled dynamic equations. 

In the dimensionless variables introduced by (|^), the dynamic equations with 
account of nonlinear terms of the second order, look: 

% + V % + <tn» + %- Re' 1 An iT 1 (0 + £g) = Cp 2 (27) 
fr + ^ff + ^-Re- 1 Aw-R- 1 (£§- z +^)=V3 

with a vector of the second-order nonlinear terms ip in the right-hand side (a 
non-dimensional value /?* = p/po is used in the right-hand side, the asterisk will 
be omitted later) : 

-"i-<4l+pf (28) 

One could rewrite a system ( p7| ) in another way: 

d 

—ip + Ltp = (p, (29) 

where the vector state (p is a specific perturbations column (vector of the fluid 
state) 









M 




i-( 




V ^3 y 






(30) 



and L is a matrix operator 



Vod/dx e 2 d/dx t 2 d/dz 

d/dx V d/dx~ Re^A- R- 1 d 2 /dx 2 cf) - R~ x d 2 jdxdz 

d/dz -R~ 1 d 2 /dxdz V Q d/dx - Re' 1 A - R- 1 d 2 /dz 2 

(31) 

All the projectors do commute with operators d/dt ■ I and L, so one can act 
by projectors on the system of equations (|30|), ( |3l| ) directly thus obtaining the 
evolution equation for the correspondent mode. There are three such equations 
for the specific perturbations p, u, w for every mode (that are equivalent) ac- 
counting the relations between these specific perturbations. So, as independent 
variable for every mode a single specific perturbation, such as pressure or one 
velocity component, may be chosen. 

Due to the existing tradition it is convenient to use the stream function as 
a basic for the TS mode and pressure perturbations for the acoustic modes. 

Acting by the Pts on the both sides of (|3C|), one gets an evolution equation: 

Ad^/dt + VoAdV/dx - d^/dx ■ d^/dz - Re^A 2 ^ = JL[-udu/dx 
—wdu/dz + pdp/dx] — -3- [—udw/dx — wdw/dz + pdp/dz] . 
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Next it should be noted that in the right-hand nonlinear side p, u, w are 
overall perturbations to be presented as a sum of specific perturbations of all 
modes: 

u = d^/dz + eA-^dpAi/dx - eA- 1 / 2 dp A2 /dx 
w = -dV/dx + eA-^dpAi/dz - eA- 1 / 2 dp A2 /dz (33) 
p = 2d 2 /dx 2 A-\^) + pax + PA2 

Indeed, there is also a density perturbation in the right-hand nonlinear side 
that was not involved in the left-hand linear one at all. The continuity equation 
reads 

dp/dt + Vodp/dx + (du/dx + dw/dz) = —d(pu)/dx — d(pw)/dz (34) 

Comparing linear left-hand side of equation (Q) with that of the first equa- 
tion from ( p6| ) , the obvious relations for the both acoustic modes follow: 
PAi = e 2 PAi , PA2 = e 2 PA2 ■ 

A limit e = yields in the TS mode for incompressible flow: pxs = . 
Therefore, the last relation for the overall density perturbation looks 

p = e 2 p A i + e 2 p A 2 (35) 

Finally, (|3^) goes to 

A* t + V Q A^ X - V X (I) Z - Re^A 2 ^ = -^ Z A^ X + * X A*,- 

e (A*A 1 /2( :PA1 - PA2 ) + Mi^-^ipAx - P A2) X + (36) 
A* z A- 1 / 2 (pAi- P A2)z) + 0(e 2 ). 

Derivatives are marked with lower indices. The first two nonlinear terms in the 
right-hand side of (^) expresses the TS mode self-action, and the last ones - 
cross acoustic-vorticity terms responsible for the acoustic mode influence on the 
TS mode propagation. The structure of the quadratic nonlinear column ( |30| ) 
yields in the absence of quadratic acoustic terms in (|36"|). 

In the limit of Vq — , <fr = the only self-action give the well-known 
evolution equation for vorticites transition jlTj] follows from (|36"|): 

A* f - Re^AH + * 2 A^ - ^ X A^ Z = 0. 

Let an acoustic field consists of only the first mode. Acting by projector Pai 
on the system (^8|), we get an evolution equation for this mode: 

dpAi/dt + V dpAi/dx + A 1 / 2 p A i/e — \[—udp/dx — wdp/dz+ 
(du/dx + dw/dz) (Zp + Spe- 2 )] + 
(-d 2 /dx 2 A- 1 <j>d/dzA- 1 + d/dx{l/2e)A- 1 / 2 ) [-udu/dx - wdu/dz + pdp/dx] + 
(d 3 /dx 3 A- 1 <f>A- 1 + {l/2e)A- 1 l 2 d/dz) [-udw/dx - wdw/dz + pdp/dz] . 

(37) 
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Here constants Z and S are denned earlier by (g). The variables p, u, w, are 
overall perturbations accordingly to (|33|), with pA2 = . So, ( pS7| ) goes to the 
final version for the directed acoustic mode 



e (d PA1 /dt + V d PA1 /dx) + A^pa, = 
i\-^zPAix + ^xPAiz] + eA- 1 ^ 2 [^ z A 1 ^p Alx + 
VxAVtpAu - 2^ xz A- 1 / 2 PAlxx - 2^ zz A' 1 ' 2 p Alxz + 2^> xx A- 1 / 2 PAlxz + 
2^ xz A- 1 ' 2 p Alzz \ + A-Va [- (^ xz f + 1^ 22 )] + 

e [-^A- 1 ^*)^ + ^xA- 1 ^*)^] + 
eA~ 1 (f>A~ 1 d 2 /dx 2 [V Z AV X - V X AV Z ] + 0(e 2 ) 

(38) 

Between the nonlinear terms one can recognize interaction (Al — TS )and gen- 
eration ones (TS — TS). The equation for p A i is obtained by projecting P A 2 and 
looks very similar. The complete system includes this equation and d3S|) , (|32|) . 
The system covers all possible processes description up to quadratic terms ap- 
proximation. Here a multimode TS waves in the (OS) equation solutions basis 
could be incorporated . The long- wave limit of such disturbance leads to coupled 
KdV system @. 

6 Resonance interaction of acoustic and TS modes. 



Equations (|36|) , ( p8| ) form a coupled system of evolution equations for interacting 
acoustic and TS modes. In the case of the TS mode generation by an incoming 
first acoustic mode, the early stage of evolution (for small amplitudes of TS 
mode) is defined by a system: 



A* t + V Q A^ X - ^ x -4> z - Re^A 2 ^ = 
-e (AM/ • A^pax + A^ x ■ A^I 2 p AVx + A* z ■ A^pau 



(39) 



e (dpAi/dt + VodpAi/dx) + A 1 / 2 p A i = § [-^zPau + ^xPAXz] + 

eA-V 2 [-H> z AV 2 P Au + * x AV 2 p AU - 2* xz A-V 2 p Alxx - (40) 
2^> zz A' 1 / 2 pai xz + 2^> xx A- 1 / 2 pai xz + 2^ xz A- 1 ' 2 pai zz \] 

All quadratic terms relating to TS — TS interaction, are left out of account . 
For simplicity, we consider only the first incoming acoustic mode. 

As it follows from the discussion in the introduction, let us find a solution 
in the form: 

Pai(x, z, t) = A\(fix, p,t)iri exp(£(cj]i — k\x)) + ^(/irr, fit)^ exp(i(ui2t — k%x)) + c.c. 

(41) 



M>(x, z, t) = B^(fix, nt)ipz(z) exp(i(w3t — k$x)) + B±([ix, fd)ipi(z) exp(z(w4i — k^x)) + c.c. 

(42) 
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where III = 7Ti(fci,u/i, z) exp(i(u>it — kix)) , II2 = ^(fe, ^2, z) exp(i(w2^ — k 2 x)) 
are planar waves. Ill satisfies the linear evolution equation ( |26| ) 

<9IIi/dt + Vo^ni/Sx + e^A^Ux = 0. (43) 

that leads to the equivalent equation for 111 

iwiHi - ikiVoIli + e _1 A 1/2 IIi = 0. (44) 

Suppose the vertical gradients of all wave functions inside the viscous layer are 
much bigger then horizontal ones. From experiments (e.g. [0), it is known, 
that wavelength of TS mode is much greater than a thickness of the boundary 
viscous layer for common values of Reynolds number. So, the operator A 1 / 2 
may be evaluated as the generalized operator (Taylor) series with respect to 
d z /k\. Hence the operator radical in the first approximation is evaluated via 
Gataux derivative as 



A^nx = ^J-kf + dlUi = ik 1 ^Ji-dl/kju 1 ^ik 1 (i-d 2 z /2kl)n 1 

and, for (^), one arrives to the ordinary differential equation, that we can 
consider as a spectral problem with the spectral parameter ki, 



(1 - ^Vo)^ + -^(1 - -L^TT! = 0. (45) 



The same equation obviously define TT2 , it is enough to change indices 1 — » 2 in 
the operator. The functions ip3( z ) > tpii 2 ) are solutions of the OS equation ( p"7| ) 
suitable for a concrete problem. A\ B4 are slowly varying functions of x, t, 
that's why an additional small parameter \i is introduced , generally, they are 
complex functions. Calculating the right-hand nonlinear expressions, we take 
only first term in series to avoid small terms of the higher order. 

Let us discuss a possibility of four- waves resonance. Examining the algebraic 
relations between parameters yields the appropriate conditions: 

UJl = 0J 2 - UI3, IV 2 = OJl + Ld 3 , CJ3 = U)\ + OJ4, UJ4 = ld 3 - LOl (46) 



Substituting the formulas Q41j , |42|) to fl43| , p9D , and picking up the resonant 
terms only, one goes to the further system of equations (complex conjugate 
values marked with asterisks, k\ — k2 + k 3 — Afc, k 3 — k\ — k^ = Afc'): 

[i ^eAiTTTi + Aix ^eVo7ri - ik\ j nidz 

eA 2 B% (1.5 (ik^* 3 ir 2z + ifeV&Ta)) 



M ^€A 2 t^2 + A 2 x (eVoW2 - ik2 J i^idz 
eAiB 3 (1.5 {-ik^l^u + ikiip az wi)) 4 
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z 



o 



A« (Baripazz + B 3X [(V + AiksRe- 1 ) t(j 3zz + (2k 3 uj 3 - 4> z ) fa]) 







(j, (B iT ip4, zz + B iX [(V Q + Aik^Re x ) ip 4zz + (2fc 4 w 4 - <f> z ) fa] ) 




We take into account that acoustic wavenumbers are real, but wavenumbers 
of both T-S modes may be complex in the general case, namely the points of 
real values form the neutral curve JlTj . 

This resulting equation may be considered as 4- wave resonance equation but 
without synchronism condition p3| . The coefficients of the equation depend on 
z, it is due to our choice of the only one transverse mode for each "horizontal" 
one. Following the lines of that have a connection of Galerkin numerical 
method, 

We continue the projecting procedure considering the transverse modes as 
a basis. In such problems we naturally arrive at two bases. One arises from 
TSW theory and its origin is from OS equation (denoted by ). Other is from 
sound problem. Of course, such bases are not orthogonal. Hence we could 
only multiply each equation by the its own basic vector and integrate across 
the boundary layer. The result is written in " back" -re-scaled variables: we put 



where the group velocities and nonlinear constants are expressed via the inte- 



Mt + c al A lx 
Mt + c a2 A 2X 
B 3 t + cts\B 3 x - 

B 4 t + CTS2B4X = 



= n al A 2 B^ Ak ^ x , 
= n a2 A 1 B*e^ Ak ^ x , 
= -eryrsiAiB 4 eW>*, 

m A* II „i( — A.k f )x 



(47) 
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grals across the boundary layer with a width 5: 



rial 



s 

f [7T? — e~ 1 iki7Ti f TT\dz\dz' 
J Voirfdz 

J^[l-5(—ik^^Tr 2 2+ik2^l x Tr2t) + f TT 1 (-ik^tp^ z Tr2 z +ik 2 tl)^ zz Tr2)dz]dz' 



Jo Vgirfdz 

5 

J [tj— £~ 1 ik 2 Tr2 f TT2dz]dz' 



C a 2 — 71 



fo Vo^jdz 



fg 7T2[(1.5(— ik3lf!3Tri z +iklJl>3 z TTl) + J (-ik 3 1p3zTr l z+ik l 1p 3z zTr l )dz)]dz' 

Ka2 = WV^, ( 48 ) 

So Vonfdz 

_ fo ■4'3[(VQ+iAk3Re- 1 )lp3zz + {'2k3U3+<Pz)lp3])dz 



CTS1 — r « , , ■ 

fo , <p3(ip4 zz ni z -k 1 k 4 ip 4zz f Tr 1 dz+ip 4zzz ir 1 )dz' 

n TSl = rS , , ^ 

Jo -4>3lp3 z zdz 

_ f 6 ijj 4 [{Va+iAk 4 Re- 1 )ip 4zz + {2k 4 u 4 + Vz )ip 4 \)dz 
° TS2 ~ f s ^ 4zz dz 

So ^ii^3 z z^l z - k l k 3i>3 zz S 7r l dz +^3 zzz 7r l) dz> 



nTS2 



S lp 4 1p 4zz dz 



A structure of the obtained equations is the particular case of general N-wave 
system , that may be solved by special technics valid for integrable equations 
p3j. The 4- wave approximation may give rise to such solutions that exhibit 
effective energy exchange between modes. The form of the nonlinearity is typical 
for a 3-wave systems and even small fluctuations of a TS field could initiate a 
rapid growth of both components if the acoustic field is big enough. It is pure 
nonlinear instability that may be supported by linear stability curve shift 
The numerical evaluation of the integrals in constants ^ need rather compulsory 
calculations. The solution of the system (^) pose also the separate problem. 
We plan to present the results in a next paper. 



7 Conclusion 

The resulting system (|3^),([38|) and the equation of the opposite directed A-modc 
could be considered as a basic one for all-perturbations over a BL description. 
We also would note that the boundary layer width may depend on x. A slow 
dependence is usually accepted and do not change the general structure of the 
expressions. As the resonant as a non-resonant processes may be studied with 
the acoustic waves separated. The modes have the different scales, hence a 
numerical modelling of the mutual generation and control also could be more 
effective. 
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